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Universal scaling law in long gamma-ray bursts 
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Abstract 

Overwhelming diversity of long gamma-ray bursts (LGRBs), discovered after the launch of Swift satel- 
lite, is a major obstacle to LGRB studies. Recently, it is shown that the prompt emission of LGRBs 
is classified into three subclasses: Type I, Type II LGRBs populating separate fundamental planes in a 
3D space defined by the peak luminosity, the duration, and the spectral peak energy, and outliers not 
belonging to either of the planes. Here we show that Type I LGRBs (LGRBs I) exhibit different shapes 
of light curves from Type II LGRBs (LGRBs II). Furthermore, we demonstrate that this classification 
has uncovered a new scaling law in the light curves of LGRBs II spanning 8 orders of magnitude from 
the prompt to late X-ray afterglow emission. The scaled light curve has four distinct phases. The first 
phase has a characteristic time scale while the subsequent three phases exhibit power law behaviors with 
different exponents. We discuss its possible interpretation in terms of the emission from an optically thick 
fireball propagating in the cricumstellar matter at relativistic speeds and argue that the observed four 
phases correspond to its hydrodynamical phases. Our classification scheme can pin down the intrinsic 
luminosities of LGRBs II through the scaling law from a sample of polymorphic GRBs. Further refinement 
of this scheme and scaling law will make a subclass of LGRBs a standard candle as reliable and accurate 
as Type la supernovae in the more distant universe than supcrnovae can reach. 
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1. Introduction 

Gamma-ray burst (GRB) is the most energetic explo- 
sion in the universe. A GRB of which the prompt emission 
lasts longer than 2 seconds is classified into a long GRB 
(LGRB). The prompt gamma-ray emission is followed by 
afterglow emission at longer wave-lengths. LGRBs are 
thought to come from the death of massive stars. The dis- 
covery of supernova 2003dh associated with a long GRB 
030329(Stanek et al. 2003) was thought to conform such 
a hypothesis. However, observational success after the 
launch of Swift satellite (Gchrels et al. 2004) casted doubt 
on the bimodal classification based on the duration of the 
prompt emission. In spite of its extremely long duration 
(~ 102 sec), GRB 060614(Gehrels et al. 2006) could result 
from a neutron star merger because of no associated super- 
nova. Even if associated supernovae are observed, many 
of them, e.g. GRBs 980425 and 060218, are far dimmer 
than ordinary GRBs. In addition to these mysteries, di- 
versity of the early X-ray afterglow discovered with the 
X-ray Telescope (XRT) on board Swift satellite is another 
major obstacle in the field of GRB study. Though all of 
these observations indicate a limitation of short-long di- 
chotomy of gamma-ray bursts, until now there has been 
no successful classification scheme to take into account all 
of these observational properties consistently. 

Well known correlations between observed quantities in 
prompt emission such as the spectral peak energy (E p ) 
- isotropic energy (Ei so ), E p - peak luminosity (L p ), 
or spectral lag (£i ag ) _ peak luminosity (L p ) correla- 
tions(Amati et al. 2002; Yonetoku et al. 2004; Norris et al. 



2000) are sometimes used as discriminators of GRBs(Lii 
et al. 2010; Gehrels et al. 2006). These correlations that 
divide GRBs into short GRBs (SGRBs), LGRBs, and 
some low-luminosity GRBs, have left a large dispersion in 
the correlation of each class. Though the dispersion had 
been implicitly assumed to follow a single gaussian dis- 
tribution in many studies before Swift and Fermi, recent 
observations of no associated supernova with an apparent 
LGRB and diversity in X-ray afterglow might have indi- 
cated that this is not the case. Recent observations, by 
Swift and Fermi etc, provide so plenty of data that this 
implicit assumption can be tested. 

Thanks to thorough studies to reduce systematic errors 
in E p , L p , Ei so and the duration of GRBs(Kancko et al. 
2006; Yonetoku et al. 2010; Kocevski 2012), criteria for the 
gold sample of GRBs were established as follows (Tsutsui 
et al. 2011; Tsutsui et al. 2012): 

1. Time integrated spectra must be fitted with the 
Band model to reduce systematic errors due to dif- 
ferent choices of spectral models. 

2. Light curves must have a time resolution of 64 msec 
to derive the peak luminosities with the same time 
resolution in the GRB rest frame by re-binning the 
fluxes. 

3. Peak photon counts must be 10 times larger than 
the background fluctuation for secure estimations of 
the duration and E[ so . 

The gold sample compiled by Tsutsui et al. (2012) had 
the smallest systematic error and made it possible to find 
out a generalized correlation between E p , TL (= Ei SO /L p ) 
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and L p . Furthermore, the correlation succeeded in dis- 
criminating, at least, three subclasses in LGRBs(Tsutsui 
ct al. 2011; Tsutsui et al. 2012) : Type I, Type II, and out- 
liers. This new classification was verified from the tempo- 
ral properties of the prompt emission (Tsutsui et al. 2012). 
Figure la shows the E p ~-T L -L p diagram of the gold sam- 
ple with the classification in Tsutsui et al. (2012). Blue 
squares indicate LGRBs I populating a fundamental plane 
plotted with the solid line, red circles LGRBs II on another 
slightly misaligned fundamental plane, and black squares 
outliers. LGRBs I and II arc also discriminated by the 
shape of their prompt emission light curves. Figure lb 
shows typical light curves for LGRBs I (990506) and II 
(081222). As exemplified in this figure, LGRBs I have 
many spiky pulses with sharp declines, while LGRBs II 
have consecutive pulses with short intervals followed by 
long tails. To clarify the difference, we make a plot of nor- 
malized cumulative counts as functions of time normalized 
with Tg§ (the period during which 98% of the total flux is 
received) for 17 LGRBs I and 11 LGRBs II in Figure lc. 
LGRBs II (red curves) deviate from the constant count 
rate (black solid line) further than LGRBs I(blue curves) 
(See Tsutsui et al. (2012) for more details). 

In this letter, we compile combined BAT - XRT light 
curves for each subclass in the gold sample and discuss dif- 
ferences in the shapes of light curves of GRBs belonging to 
different subclasses in section 2. We analyze light curves 
of LGRBs II in particular to seek universal features cover- 
ing all the phases from the prompt to late X-ray afterglow 
emission in section 3. Finally, we present a possible inter- 
pretation using a fireball model in section 4. 

2. X-ray afterglow 

The classification scheme was developed with informa- 
tion only of the prompt emission. Nonetheless, it prob- 
ably has a link with the afterglows. Here we investigate 
the temporal profile of combined prompt - X-ray afterglow 
emission for each subclass of LGRBs. To do this, we ob- 
tain combined BAT- XRT flux light curves in the energy 
range of 15-50 kcV for our gold sample from the Swift 
Burst Analyser website 1 (Evans et al. 2010). We choose 
BAT data binned with the signal-to-noise ratio of at least 
5. Swift/XRT observed 28 out of 36 LGRBs in the gold 
sample. The present sample with Swift observations is 
composed of 10 LGRBs I, 10 LGRBs II, and 8 outliers. 
In Figure 2, we show combined BAT-XRT light curves of 
all the 28 LGRBs. The reddish points indicate the light 
curves of LGRBs II, the blueish points LGRBs I, and the 
blackish points outliers. Roughly speaking, the data for 
the first 100 sec are taken with the BAT and the rest 
with the XRT. As the figure shows, LGRBs II have three 
breaks in their light curves in common, while LGRBs I 
and outliers exhibit much complicated behaviors. We will 
discuss the early X-ray afterglow emission of LGRBs I 
and outliers in separate papers. Here we concentrate on 
the analysis of LGRBs II, which is much simpler than the 



other types. 

3. Scaling law 

To seek a universal scaling relation in the light curves 
of LGRBs II, we normalize 15-50 keV luminosity Lx with 
L p and the time T = T b s /(1 + z) measured from the BAT 
trigger in the GRB rest frame with Tl, where T b s denotes 
the observed time, z the redshift of the GRB host galaxy. 
In Table 1, we summarize redshifts, classes, Xl, and L p 
for the gold sample observed by Swift/XRT. In Figure 3a, 
we make a plot of normalized luminosity Lx/L p versus 
normalized time T/Tl for LGRBs II. We find that two 
of them, showed in the inset, have strong X-ray flares 
and exclude them from the following analysis. These two 
events also seem to follow the same scaling law for the 
others, except for X-ray flares. This suggests that the 
existence of X-ray flares does not affect the rest of the 
afterglow emission. 

We found that the lack of data between the first and 
second breaks for some GRBs makes it difficult to fit a 
function to the light curve of each LGRB II. Instead of 
fitting to data of each LGRB II, we stack data of the 8 
LGRBs II thus selected to investigate the average proper- 
ties and fit a broken power law model with three breaks. 
We use the least square method without any weighting 
factor, because the dispersion of the data around the best 
fit model is larger than error bars. Therefore we minimize 
the following least square merit function, 
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and define the standard deviation as a = y/S m i n /(N — 7). 
This model has 7 parameters : the normalized luminosity 
^norm m the prompt (first) phase, the power law exponent 
for the i-th phase, on (i = 2, 3, 4), the first break time, tbi, 
second break time, ib2 and the third break time ib3- As a 
result, we obtain the best fit model, 
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where C x = 0.173, C 2 = (15.8/0.48)- 1 " , and C 3 = 
(273.6/15.8)- 049 . The best fit values (1-cr errors) for these 
parameters are L norm = 0.173(0.0058), a 2 = -1.99(0.026), 
a 3 = -0.49(0.038), a 4 = -1.31(0.018), t hl = 0.48(0.013), 
t h2 = 15.83(0.92), and t h3 = 273.6(27.6) with standard de- 
viation a = 0.84. We find that the fourth phase has a 
slightly larger standard deviation (174 = 1.14) than the 
other phases (ci = 0.68, cr 2 = 0.68, and CT3 = 0.65) where 
Oi is the standard deviation for the i-th phase. Jet breaks 
might be responsible for this large dispersion in the fourth 
phase. In Figure 3a, we show the best fit model with the 
solid red lines. 
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Fig. 1. Classification scheme of prompt emission of long GRB. a, The E p -Tl-L p diagram for 36 long GRBs in the gold 
sample in Tsutsui ct al. (2012) with the residual plot in the bottom panel. Three subclasses based on the distribution in this plot 
are indicated by different colors: Type I (blue), Type II (red) and outliers (black). The same color is used for the same subclass 
in the other panels. The solid black line indicates the fundamental plane for Type I long GRBs. b, Examples of normalized light 
curves, c, Normalized cumulative light curves. The black solid line indicates a constant count rate. 



Table 1. Spectroscopic rcdshifts, class, Tl, L Pj 5i(=L p /(10 51 erg/s)) of the gold sample observed by Swift/XRT from Tsutsui et al. 
(2012). To minimize the dispersion in the Ep-T^-Lp diagram, L p were re-binned by 2.752 sec in the GRB rest frame. 



GRB 


z 


class 


[sec] 


Lp,51 

[erg/sec -1 ] 


GRB 


z 


class 


Tl 

[sec] 


Lp,51 

[erg/sec -1 ] 


050401 


2.9 


Type I 


7.00 


7.34 


090424 


0.544 


Type II 


3.97 


4.27 


050525 


0.606 


Type I 


7.95 


71.80 


090618 


0.54 


Type I 


21.7 


69.05 


050603 


2.821 


Type II 


3.51 


227.31 


090902B 


1.822 


Type I 


6.37 


5.82 


061007 


1.261 


Type I 


9.13 


81.96 


090926A 


2.106 


outlier 


5.00 


334.17 


070125 


1.547 


Type I 


8.60 


33.29 


091003 


0.897 


outlier 


8.40 


62.85 


071003 


1.604 


outlier 


8.35 


17.10 


091127 


0.49 


Type II 


8.45 


67.297 


071010B 


0.947 


Type II 


5.50 


219.39 


091208B 


1.063 


Type I 


4.41 


299.80 


080319B 


0.937 


outlier 


19.26 


1.16 


100414A 


1.368 


outlier 


13.0 


7.52 


080413B 


1.1 


Type II 


3.18 


44.17 


100906A 


1.727 


Type II 


9.13 


11.773 


080721 


2.602 


Type I 


3.59 


2.92 


110213A 


1.46 


outlier 


4.80 


11.9 


081121 


2.512 


Type I 


3.95 


135.53 


110422A 


1.77 


outlier 


4.54 


483.4 


081222 


2.77 


Type II 


3.14 


107.56 


110503A 


1.613 


Type II 


2.96 


384.026 


090323 


3.57 


Type I 


12.84 


107.14 


110715A 


0.82 


Type II 


4.42 


10.905 


090328 


0.736 


outlier 


22.61 


48.96 


110731A 


2.83 


Type II 


2.39 


2.0 
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Fig. 2. Combined BAT-XRT light curves for each subclass of the gold sample, a, Overall shape of the light curves of 
GRBs belonging to each subclass (blueish, Type I: redish, Type II: blackish, outliers). Individual events in the same subclass are 
distinguished by gradation, b, Close-up of the region enclosed by a rectangle in the left panel. 
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Fig. 3. Normalized (a) and unnormalized (b) BAT-XRT light curves for Type II long GRBs. Luminosities normalized 
with L p as functions of the normalized time T/Tl are shown in the panel a. The panel b shows light curves constructed from the 
same data as in the panel a. Two events with strong X-ray flares are separately plotted in the inset. The best fit curve is indicated 
by the red solid line with three phases described by power laws with different exponents (see text for details). 
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Table 2. Relationship between time scales, brightnesses and 
standard deviations for different normalizations. The normal- 
ization adopted in equation (2) is separately listed in the last 
row. 



Time 


Luminosity 


a 


-* obs 

T 

T/T 90 


Fx 

Lx 

Lx/L p 


0.93 
1.02 
1.01 


T/T L 


Lx/L p 


0.84 



To demonstrate how our scaling law works well, we fit 
the same model to three light curve plots constructed 
from the same data with different normalizations, i.e., 
Fx(T bs), Lx{T), and L x (T ohs /T 90 )/L p . The standard 
deviation in each fitting in table 2 shows that the normal- 
ization in equation (1), i.e., Lx(T bs/TL) / L p reduces the 
dispersion. We summarize the standard deviation in each 
fitting in table 2. As table 2 shows, the normalization in 
equation (2), Lx/L p versus T/Tl, improves the fit of the 
model. For a visual comparison, we show a relationship 
between the intrinsic luminosity Lx and time T in Figure 
3b. Figures 3a and 3b show that our normalization really 
reduces the dispersion of light curves. 

4. Implication &: discussion 

Here we briefly discuss a possible interpretation of the 
universal scaling law presented above by using the tempo- 
ral evolution of photosphcric emission of an optically thick 
jet propagating in the circumstcllar matter (CSM). The 
dynamical evolution of such a jet relevant to observations 
is composed of four phases: the jet emergence phase, the 
stationary expansion phase, the blast wave phase, and 
the reverse shock phase. The last two phases are de- 
scribed by self-similar solutions (Nakamura & Shigeyama 
2006; Shigeyama et al. 2012). In the first two phases, 
the central engine must be still active and the duration 
of the activity determines the time scale Tl. The pre- 
dicted temporal evolution of the luminosity in each phase 
is summarized as follows. The observed light curves in the 
first two phases seem to suggest the temporal evolution of 
the X-ray luminosity in the form of Lx oc (1 +T/Tl) -2 . 
Furthermore, the exponents in the other phases can be 
deduced from the corresponding self-similar solutions for 
propagation of relativistic shocks in a stationary wind: 
Q!3 = 1 — V3, «4 < —3/2. Detailed derivation is given in 
a separate paper (T.S. and R.T. in preparation). The 
rough agreement between these exponents and those in 
equations (2) deduced from observations encourages such 
a modeling. 

Deeper understanding of the scaling law and refinement 
of the classification scheme will lead to more accurate de- 
termination of the intrinsic luminosities of LGRBs and 
then shed light on the more distant universe than ever 
reached by supernovae. 
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